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Abstract. The nonlinear kinetic model of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used to 
describe the relevant properties of Cassiopeia A (Cas A). In order to reproduce the SNR's observed size, expansion 
rate and thermal X-r ay emission we e mploy a piecewise homogeneous model for the progenitor's circumstellar 
medium developed bv lBorkowski et al.l dl996l) . It consists of a tenuous inner wind bubble, a dense shell of swept-up 
red supergiant wind material, and a subsequent red supergiant wind region. A quite large SNR interior magnetic 
field Bd ~ 1 mG is required to give a good fit for the radio and X-ray synchrotron emission. The steep radio 
spectrum is consistent with efficient proton acceleration which produces a significant shock modification and leads 
to a steep electron spectrum at energies e c < 1 GeV. The calculated integral 7-ray flux from Cas A, F 1 tx e" 1 , is 
dominated by 7r°-decay 7-rays due to relativistic protons. It extends up to roughly 30 TeV if CR diffusion is as 
strong as the Bohm limit. At TeV energies it satisfactorily agrees with the value 5.8 x 10~ 13 cm~ 2 s -1 detected by 
the HEGRA collaboration. 
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1. Introduction 

Cassiopeia A (Cas A) is a prominent shell type su- 
pernova remnant (SNR), and a bright source of syn- 
chrotron radjat ion_obserro 

Bell et alJll973 lTuffslll98fil lBrannlll987t lAndersen et alJ 
199ll iKassim et al.lll995|) and most probably also in the 
X-rav band (|Allen et al.lll997t iFavata et alJll997^ . Since 
SNRs are widely suspected to be the main sources of 
the Galactic cosmic rays (CRs) up to energies of at least 
10 15 eV, this youngest Galactic SNR offers an important 
possibility to test this CR origin hypothesis. With the syn- 
chrotron interpretation, the measured emission of hard X- 
rays is direct evidence for the existence of a large num- 
ber of relativistic electrons with energies up to about 
10 TeV, presumably accelerated at the outer SNR shock. 
Experimental information about nuclear CR production 
can only be obtained by high energy 7-ray measurements. 
If protons are accelerated in Cas A to at least the same 
energy and as efficiently as electrons then the 7r°-decay 
7-ray spectrum, created in their hadronic collisions with 
the background nuclei, should extend to energies above 
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1 TeV with a hard power-law. The detection of a signal 
in TeV 7-rays has bee n indeed recently repo rted by the 
HEGRA collaboration l|Aharonian et al.ll200l[) . Therefore 
it is of the exceptional importance to find out whether this 
7-ray emission is consistent with a hadronic origin. 

In this paper a thorough th eoretical analysis wi ll 
be given, based on ear lier work feerezhko et alJ 1200 lib 
It fo llows similar analyses f or SN 1006 llBerezhko et alJ 
|2002|) and for Tycho's SNR l)Volk et al.ll2002jb We note 
that in the meantime also new 7-ray observations of 
SNR RXJ1713. 7-394 6 by the CANGARO O experiment 
have been discussed i Enomoto et alJl2002ti and criticized 
l|Reimer fc Pohlll2002HButt et alJl2002h on a phcnomcno- 
logical basis. Unfortunately, the latter object is up to now 
only very poorly understood, for example regarding the 
type of explosion (quasi-explosive nuclear burning of an 
accreting White Dwarf vs. core collapse of a massive pro- 
genitor star) and the distance (estimates range from 1 kpc 
to 6 kpc). The object also lies in an environment of com- 
plex morphology. In addition the radio observations are 
so scarce that they do not allow to determine a spec- 
trum. As a result only phenomenological single box-type 
estimates for the particle spectra which are responsible 
for the synchrotron/Inverse Compton radiation and the 
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hadronic 7r°-decay 7-ray emission have been made up to 
now, preventing any firm conclusion regarding the nature 
of the 7-ray emission. Even a detailed theoretical model 
like ours could, with the currently available data, only lead 
to a number of constraints. In contrast, CasA is perhaps 
the best-studied object of its kind. The following analy- 
sis makes detailed use of this impressive multi- wavelength 
knowledge. 

The chemical compositions of characteristic parts 
of t he supernova (SN ) ejecta (the fast moving knots, 
e.g. iReed et alJ Il995l and references therein) and of 

t he circumstell ar ma terial (the quasi-stationary floc- 

culi. IPeimbertl Il97lt I Pcimbcrt fe van den Berghl Il97ll 



iKirshner fe Chevalier Il977t I Chevalier fe Kirshnerl fl97i 
suggest that the progenitor star started as a massive ob- 
ject. It should have evolved from a hot main sequence star 
with a tenuous fast wind to a red supergiant phase with 
a slow but dense wind, and finally into a hot Wolf-Rayet 
star with a fast wind again, sweeping up part of the pre- 
ceeding red supergiant material into a dense shell before 
exploding as a core collap se SN llChevalier fc Lianell 19891 
iGarcia-Seeura et alJll996 ). In this spirit iBorkowski et all 
l)l996|) modeled the thermal X-ray emission as well as 
the size and the expansion rate of the present SNR 
from their result that the outer blast wave had already 
passed the swept-up red supergiant wind shell and was 
presently propagating through the unperturbed slow red 
su pergiant wind. Acc ording to the X-ray measurements 
of iFavata et alJ l)l997j) the SNR shock has swet up about 
5Mq of this circumstellar material. Therefore we can ex- 
pect a detectable level of 7r°-decay 7-ray flux. 

Empirical arguments for such a structure of the cir- 
cumstellar material around the presupernova star come 
for example from the radio observations of CasA. Despite 
its youth of about 320 years, characterized also by the rel- 
atively small amount of matter swept up by the SN shock, 
the obser ved radio flu x undergoes considerable secular de- 
cline (e.g. iReeslllQQnfr . which amounts to about 0.6% per 
year. For a uniform circumstellar material one would ex- 
pect the radio flux to increase until an intermediate Sedov 
stage of evolution, corresponding to an age of about 10 3 yr. 
Only in the case when the SN shock propagates through a 
wind of the progenitor star with progressively decreasing 
density p cx r~ 2 , the to tal radio emission is expected to 
decrease with time (e.e. iBerezhko fc Ksenofontovl ll999). 

This suggests two rather complementary approaches 
regarding the production of nonthermal particles and their 
radiation. One can either emphasize the strong radio emis- 
sion from small inhomogeneities whose free energy (ki- 
netic and magnetic) is assumed to lead to stochastic 
acceleration of energetic electrons which are responsible 
for the s ynchrotron emission from the radio to the X- 
ravhaud llScott fc ChevalierjT97ilDickel fc Oreiserll 9791 
arkarl [ 



does not have to play an important role for electron ener- 
gization, but it may still be the ma in source for the non - 
thermal nuclear particle component I Atovan et all2000ah . 
An inventory of the radio to high energy 7-ray emission 
along these lines has b een made phenomenologically by 
lAtovan et al.l l)2000bla[) . including detailed fits to the ob- 
served synchrotron spectrum. 

Alternatively one can to lowest order ignore the role 
of small-scale inhomogeneities for the production of the 
very high energy nonthermal particle component. At least 
as its carriers they are in difficulty on the argument that 
high energy particles (ultrarelativistic electrons and nu- 
clear particles) require a large acceleration/propagation 
volume given by spatial scales L ~ A m f p (p) x c/u. Here 
u denotes a characteristic nonrelativistic speed of mass 
motions, c is the speed of light, and A m f p (p) is the scat- 
tering mean free path that increases with momentum p. 
For strong scattering A m f p = 0(r g ), where r g oc p is the 
particle gyro radius. At least for the highest energy par- 
ticles L amounts to a few percent of the remnant radius. 
Ignoring the dynamical effects underlying the small-scale 
emission features, the dominant remaining large-scale en- 
tropy generator is the SNR blast wave. It has then to be 
investigated which of the observed nonthermal features 
can be explained by shock accelerated particles. 

In this paper we shall apply the nonlinear ki- 
netic model for diffusive shock acceleration in SNRs 



Cowsik fc Sarkari 11984 jjones et all Il994t lAtovan et alJ 
2000b). This gives also rise to an attendant nonthermal 
Bremsstrahlung (NB) and Inverse Compton (IC) compo- 
nent extending to the high energy 7-rays. In this phe- 
nomenological scenario the large scale SNR blast wave 
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the specific case of CasA. This spherically symmetric 
model should describe the global properties of the non- 
thermal emission especially at high energies, even though 
small-scale processes may indeed introduce some modifi- 
cations at low particle energies, contributing especially to 
the radio flux. The calculation of the nonthermal particle 
populations is a central part of the theory which also de- 
termines the spherically symmetric portion of the overall 
gas dynamics in the remnant selfconsistently. 

In fact, we shall attempt to investigate whether the 
global properties of Cas A are consistent with the idea that 
the SN blast wave is the main source of energetic particles 
and in particular, whether the observed 7-ray emission is 
consistent with a hadronic origin. 



2. Model 

To descr ibe the circumste l lar m edium we use the specific 
model of IBorkowski et all l|l996|) . Accordingly, part of the 
slow red supergiant wind of the SN progenitor has been 
swept up into a dense shell by a fast stellar wind during the 
final blue supergiant (probably Wolf-Rayet) phase of the 
progenitor star. Therefore the inner circumstellar medium 
consists of three zones: a tenuous wind-blown bubble, a 
dense shell, and a freely expanding red supergiant wind. 
The outer circumstellar regions that were generated dur- 
ing the main sequence phase play no role here. 
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We describe the profile of gas number density N g = 
p/rrip in the analytic form 



, T N hsh + N w iVw-AW ifr-R 2 
N e = + tanh ( — — 

where 

N h + N sh N sh -N h ,(r-R x 
JVbsh = n 1 n tanh 



and 

N w = N w2 (R 2 /r) 2 



(1) 



(2) 



(3) 



is the gas number density of the free red giant wind (region 
r > R 2 ), Nb and iV s h correspond to the hot bubble (region 
r < Ri) and shell (region i?i < r < R 2 ) respectively, p is 
the gas density, m p is the proton mass. This interpolation 
formula provides a smooth transition between the above 
three zones on the scale I which is taken small enough, 
I -C that its concrete value does not influence the 
final results. 

We use the same type of formula (1) and (2) for the 
magnetic field profile Bo(r) with B^, B s h and 



B w — B w2 R 2 /i 



(4) 



for the bubble, shell and wind regions respectively. 

The SN explosion ejects an expanding amount of mat- 
ter with energy E sn and mass M e j . During an initial period 
the ejecta have a wide distribution in radial velocity v. The 
fastest part of the se ejecta is descri bed by a power law 
dMej/dv oc v 2 ' ~ k l) Jones et al.lll98l|) . The interaction of 
the ejecta with the circumstellar medium creates a strong 
shock which will diffusively accelerate particles. 

The acceleration model consists in a selfconsistent so- 
lution of the CR transport equation toget her with the gas 
dynamic equations in spherical symmetry l|Berezhko et alJ 
1996; B erezhko fc Volkl Il997|) . in an extension of this 
model to the case of a no nuniform circumstellar medium 
^Berezhko fc Volkll2000h|). 

The CR diffusion coefficient is taken at the Bohm limit 



K(p) = n(mc)(p/mc), 



where re(mc) 



— 



(5) 

/(3eB); e, m and p are the particle 
charge, mass and momentum, respectively; B is the mag- 
netic field strength; and c is the speed of light. 

In the downstream region the magnetic field is assumed 
to be frozen into the gas and is described by the equation 



dB 

~dt 



V x (w x B), 



(6) 



where w is the gas mass speed, directed everywhere radi- 
ally in our spherical model. The boundary condition for 
this equation is the relation between the postshock mag- 
netic field B 2 and the upstream field at the shock position 
B s = Bq(R s ). The tangential field component B± is am- 
plified at the shock front because of compression, so that 
its postshock value is B 2 ± = <tBq±, where a is the shock 
compression ratio. Therefore we have B 2 = ctb-Boj where 



ob = \J (c 2 — 1) sin 2 4> + 1, <fi being the angle between the 
upstream magnetic field Bq(R s ) and the shock normal. 
For simplicity we suggest that the downstream magnetic 
field is purely tangential and use ctb = cr. 

The number of suprathermal protons injected into the 
acceleration process is described by a dimensionless injec- 
tion parameter r\ which is a fixed fraction of the ISM par- 
ticles entering the shock front. For simplicity it is assumed 
that the injected particles have a velocity four times higher 
than the postshock sound speed. Unfortunately there is 
no complete selfconsistent theory of a collisionless shock 
transition, which can predict the value of the injection 
rate and its dependence on the shock parameters. For the 
case of a purely parallel shock hyb rid simulations pre- 
dict quite a high ion i njection (e.g. IScholer et all Il992t 
iBennet fc Ellison! Il995h which corresponds to the value 
rj ~ 10~ 2 of our injection parameter . Such a high injection 
is consistent with a nalytical theory l|Malkov fc V51klll995L 
ll99filMaikovl l 1998|) an d confirmed by measurem ents near 
the Earth's bow shock l|Trattner fc Scholerll99^ . We note 
however that in our spherically symmetric model these re- 
sults can only be used with some important modification. 
The circumstellar magnetic field may be assumed to be 
strongly perturbed at least in the shell region. Efficient 
particle injection takes place only on those portions of 
the shock surface which are currently locally quasiparal- 
lel. On other parts of the shock where it is more and more 
oblique, the magnetic field essentially suppresses the leak- 
age of suprathermal particles from the downstream re- 
gion back upstream ([Ellison e t all Il995t iMalkov fc V51kl 
119951) . Therefore the mean injection rate and subsequent 
CR production efficiency, properly averaged over the en- 
tire shock surface, is expected to be considerably lower 
compared with the case of a purely parallel shock. Due to 
this fact the number of accelerated CRs, calculated within 
our spherically symmetric model, has to be corrected (de- 
creased) by some renormalisation factor / rc < 1, because 
of the lack of efficient injection/acceleration on the part 
1 — /rc of the actual shock surface. 

We assume that electrons are also injected into the 
diffusive shock acceleration process still at nonrelativistic 
energies below m c c 2 . Since the e lectron injection mecha - 
nism is not very well known (e.g. IMalkov fc Drurvll20nih . 
for simplicity we consider their acceleration starting from 
the same momentum as protons. At relativistic proton en- 
ergies they have exactly the same dynamics as the protons. 
Therefore, neglecting synchrotron losses, their distribution 
function at any given time has the form 



Up) = K op f(p) 



(7) 



for relativistic proton energies, with some constant fac- 
tor K cp which is about 10~ 2 for the average CRs in the 
Galaxy. We take it here as a parameter and use relation 
(|7|) only at the injection momentum in order to deter- 
mine the number of injected/accelerated electrons for a 
given proton injection rate. At all other momenta the elec- 
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tron distribution function is found by solving the transport 
equation. 

In fact the calculated electron distribution function 
f e ip) deviates from expression J7J only at sufficiently 
large momenta because of synchrotron losses, which are 
taken into account by supplementing the ordinary diffu- 
sive transport equation with a loss term: 

(8) 



dt J 



~3 



dp 



1 d 



P_ 
n 



where the first three terms on the right hand side of this 
equation describe diffusion, convection due to the mass 
velocity w of the gas, and adiabatic effects, respectively. 
The synchrotron loss ti me in the third term is determined 
by the expression (e.g. iBerezinskii et alJll990^ 



n 



(9) 



where m e is the electron mass and r the classical electron 
radius. 

The solution of the dynamic equations at each instant 
of time yields the CR spectrum and the spatial distri- 
butions of CRs and gas. This allows us to calculate the 
expected flux F^(e 7 ) of 7-rays from 7r°-decay because 
of hadronic (p-p) co llisions of CRs with the gas nuclei. 
Followin g the work of Dcrmcr ( 1986) and its later improve- 
ment bv lNaito fc Takaharal l)l994|) we use here the isobar 
model at the proton kinetic energies < 3 GeV and the 
scaling model at > 7 GeV with a linear connection be- 
tween 3 and 7 GeV. This model agrees very well wit h the 
sim pler approach, introduced bvlDrurv et alJ J l994j) (see 
also IBerezhko fc Volklll997l l2000al iBerezhko et aljll999|) 
at high energies e 7 > 0.1 GeV, except in the cutoff region, 
where the scaling model yields a significantly smoother 
turnover of the 7-ray spectrum at somewhat lower ener- 
gies. 

The choice of K ep allows us to determine in addition 
the electron distribution function and to calculate the as- 
sociated emission (for details, see IBerezhko et al . 2002, m 
a recent analysis of SN1006). We calculate here the IC 
radiation taking into account as target photon fields - be- 
sides the cosmic microwave background - the infrared field 
with a mean photon energy of e p h = 0.01 eV and an energy 
density of 0.2eV/cm 3 , and the optical field with a mean 
photon energy o f e n h = 1.5 eV and an energy density of 
0.5eV/cm 3 fe.g. lDrurv et alJll994 fcaisser et alJll99S|) . 

3. Results and Discussion 

3.1. Supernova dynamics 

In the current approach, we have reduced the circumstcllar 
density profile everywhere by a factor of 1 . 8 com pared 
with what was derived by iBorkowski et alJ l|l996|) . This 
corresponds to the parameter values: 

N h = 0.17cm" 3 , N sh = 11cm" 3 , N w2 = 2cm" 3 , 
i?i = 3.8 x 10 18 cm, R 2 = 4.94 x 10 18 cm. 



in 
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Fig. 1. Shock (contact discontinuity) radius R s (R p ) and 
shock (contact discontinuity) speed V s (V p ) (a); total 
shock (a) and subshock (er s ) compression ratios (b); ejecta 
(E C j), CR (E c ), gas thermal (E gt ) and gas kinetic (E g k) 
energies (c) as a function of time. The vertical dotted line 
corresponds to the current evolutionary stage. The cir- 
cumstellar gas number density A r g (i? s ) and magnetic field 
B n (R s ) profiles are shown in (a) by dashed and dash- 
dotted lines, respectively. Scale values are R\ — lpc, 
V\ = 30 000km/ s, £q = 31. 7 years. The observed size 



( Reed et al.lll995h and speed l|Andersen fc Rudnicklll995]) 

of the shock are shown by the full circle and triangle, re- 
spectively, in (a). 



Together with the renormalization factor / ro which is dis- 
cussed later, the lower circumstellar density is the main 
reason for the significant reduction of the expected 7- 
ray flux (see below) compared to our previous calcula- 
tions which did not take into acco unt f rs and were base d 



a swept- up mass value of 8Mq (Berezhko et alJEoOll) . 
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The density reduction is close to that suggested by the X - 
ray analysis of the BeppoSAX data llFavata et alJ 19971). 
which yielded 5 Mm i nstead of 8M ? ijVink et al. I Il996t 
iBorkowski et al.lll996ft for the swept-up mass. 

In order to reproduce the observed shock size R s and 
its expansion rate V s the following SN parameters are 
used: explosion energy E sn = 4 x 10 50 erg, ejecta mass 
M e j = 2Mq, and powe r-law index k = 6 for the ejecta ve- 
locity distribution (cf. iBerezhko et alJl2002l). As dis tance 
to CasA we adopt d = 3.4 koc fcf. lReed etaTlll995() . 

The results of our calculations together with the ex- 
perimental data are shown in Fig. HUH In Fig.^L we also 
show the profiles of the gas number density N g (R s ) and 
magnetic field B (R S ), upstream of the shock front that is 
at position R s (t). 

As can be seen in Fig.^,, after the shock has entered 
the shell region its speed drops during the initial 60 years 
by more than a factor of ten and then remains almost 
constant up to the current epoch. 

At the current epoch t 320 yr th e calculations rea- 
sonably reproduce t he observed size llReed et alJ 1 19951) 
and expansion rate l) Andersen fc Rudnickl Il995|) of the 
remnant. 

A proton injection rate r\ = 2.5 x 10~ 3 is used in or- 
der to provide the nonlinear shock modification required 
to reproduce quite a steep radio emission spectrum (see 
below). According to Fig.^J) the shock is indeed strongly 
modified by the CR backreaction: the total shock compres- 
sion ratio cr w 5.5 exceeds the classical value 4, whereas 
the subshock compression ratio is considerably smaller, 
cr s m 2.7. 

3.2. The accelerated CR spectrum 

According to Fig.^: about 15% of the explosion energy 
has been transformed into CRs at the current stage, that 
is E c w 6 x 10 49 erg. 

The overall momentum spectrum of accelerated pro- 
tons 



N(p,t) 



16ttV 



drr 2 f(r,p, t), 



m p c and 



presented in Fig.[3 turns over at about p = 10 5 
extends almost up to p max = 10 6 m p c because of the ex- 
tremely high magnetic field strength used: £? s h = 200 /iG, 
B W 2 — 100 fiG. The value p max is limited by geometrical 
factors which are the finite size and speed of the shock, its 
deceleration an d the adiabatic c ooling effect in the down- 
stream region IjBerezhkol Il996j) . If we approximate the 
overall proton spectrum by a power-law N cx p~ 7 then the 
power-law index slowly decreases from 7 ft! 3 at pZSm p c 
to 7 w 1.9 at the highest momenta 10 2 m p ^p^ 10 5 ra p c 
due to the strong shock modification. 

The shape of the overall electron spectrum deviates 
from proportionality to the proton spectrum N c (p) cx 
N(p) at high momenta p > p\ « 30 m p c, because of the 
synchrotron losses in the downstream region where the 
magnetic field Ba « B2 is of the order of 1 mG. According 
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Fig. 2. The overall CR proton (solid line) and electron 
(dashed line) spectra as function of momentum. 

to expression @ the synchrotron losses become important 
for electron momenta greater than 



Pi 

m D c 



1.3 



10 s yr\ /10/iG 



t 



B t] 



(11) 



Substituting the SN age t = 320 yr into this expres- 
sion we have p\ « 30 m p c, in good agreement with the 
numerical results (Fig.|5J). 

The shock constantly produces the electron spectrum 
f c cx p~ q with q fa 4 up to the maximum momentum 
-Pmaxj which is much larger than p\. Therefore, within the 
momentum range p\ to p„ lax due to the synchrotron losses, 
the electron spectrum is f c cx p~ 5 and the corresponding 
overall electron spectrum is iV e cx p~ 3 . 

The maximum electron momentum can be estimated 
by equating the synchrotron loss time @ and the accel- 
eration time 



(10) r a = A(- + - 
Au 



«i 
til 



1l 

U2 



(12) 



where u\ = V s and u-i = u\/a are the upstream and down- 
stream gas velocities with respect to the shock front, and 
Am = u\ — U2- In the general case the downstream mag- 
netic field equals B2 — <JbBq and correspondingly 



6.7 x 10 4 



Vs 



10 3 km/s 



' (^-1) 
a(l + o- b <j) 



10 [iG 
B 



(13) 



The main part of the electrons with the highest energies 
e<J 10 TeV is produced during shock propagation through 
the shell. At this stage V s ~ 2500 km/s which leads to a 
maximum electron momentum p° lax w 10 m p c (for sim- 
plicity we use (Tb = cr below), in agreement with the nu- 
merical results (Fig.0. 
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According to the results from the next section, the 
energy content of the electrons at the current epoch is 



K cp E c 



2.4 x 10 47 erg. 



3.3. The synchrotron emission from CasA 

The parameter K cp = 4x 10~ 3 gives reasonable agreement 
between the calculated and the measured synchrotron 
emission in the radio- and X-ray ranges, as one can see 
from Fig.0 The calculated synchrotron spectrum at the 
epoch 1970 (solid curve) , which correspond s to that of the 
radio measurements of iBaars et al.l l|l977|) . is compared 
with the experimental data. Also a calculated curve for 
the present epoch (dashed line), corresponding to the X- 
ray measurements, and a prediction for the epoch 2022 
(dashed-dotted line) are given. The strong shock modifi- 
cation leads to a steep spectrum N c oc p~ 2 6 of sub-GeV 
accelerated particles (Fig.|2J. As result of our choice of a 
large magnetic field (_B s h = 200 //G, B w2 = 100 /iG) the 
calculation fits the radio data S v oc v~ a , a ~ 0.8, rather 
well. In fact, the electron spectrum has a concave shape 
at p < p\ w 30 m p c. This leads to a flattening of the syn- 
chrotron spectrum S„(v) at v = 10 to 100 GHz, consistent 
with the experiment. Our calculated sync hrotron flux fits 
even marginally the iMezger et alJ ||198(tI) measurements 
at 1.2 mm and th e infrared emission measured at 6 /^m by 
iTuffs et all l)l997ft , even though, as one can see from Fig.|3| 
these two points are above the pure power-law extrapola- 
tion S v oc i/~ - 77 . It is possible that the far i nfrared energ y 
flux has a significa nt t hermal comp onent, cf. lBraur] l|l987h . 
ITuffs et all l|l997l) and lVinkl l|l999h . The steepening of the 
electron spectrum at p > p\ caused by synchrotron losses 
leads to a flat spectral energy distribution uS u that con- 
nects the nonthermal mid-infrared with the X-ray band 
(Fig-EJ. 

It is important to note that the complicated shape of 
the electron spectrum is naturally produced by the model. 
As far as the overall nonthermal spectrum is concerned, 
there is no strong need to add other sources of energetic 
particles. The bright 'radio ring' is in any case explain- 
able by the swept-up red supergiant wind shell which is 
contained in our model. It is another question whether 
individual 'radio knots' are local acceleration regions or 
rather just magnetic field enhancements bathed in a per- 
vasive particle background. For the nonthermal X-ray and 
TeV 7-ray emission they appear too small to locally ac- 
celerate the radiating multi-TeV particles in a stochastic 
manner. For a given frequency, a radio knot that con- 
sists of a local B-field enhancement will be illuminated by 
lower energy electrons than the ambient spatial regions. 
The strong nonlinear modification of the SN shock that 
generates such particles then ensures that these lower en- 
ergy electrons have a steeper spectrum - just what distin- 
guishes these knots form the average emission. Whether 
this effect is quantitatively sufficient, is a difficult ques- 
tion on the characteristics of the turbulent structure of the 
shocked red supergiant wind shell. However, it might go a 
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Fig. 3. Overall synchrotron spectral energy distribu- 
tions as a function of frequency, calculated for the 
epochs 1970 (solid curve), 2002 (dashed curve), and 2022 
(dashed-dotted cur ve). The radio-emission above 100 MHz 
llBaars et al.lll977ft. the data at 1.2 mm (tri angle) from 
IMezger et all l|l98q) and 6/Ltm (square) from lTuffs et al " 
(I1997T) . as well as the hard X-ray spectrum IjAllen et al 
Il997j) are presented. Due to the non-contemporaneous ob- 
servations, the radio data should correspond to the solid 
curve, whereas the X-ray data must be compared with the 
dashed curve. In the radio range theory and experiment 
agree fairly well. To achieve agreement also in the X-ray 
regime, a small increase of the overall electron cutoff mo- 
mentum p ma x of roughly 10% would be required (see text 
for details). 



long way to explain the 20% contribution l|Tuffslll986|) of 
the radio knots to the total radio emission. Potential 2 nd 
order Fer mi acceleration in the t urbulent shell, as pro- 
posed bv IScott fc Chevalierl l|l975|) . cannot be ruled out. 
However, there is no indication for its action either. For 
diffusive shock acceleration at the blast wave, on the other 
hand, we know from theory that it can have the required 
high efficiency. And it can do it all. No other process needs 
to be invoked. 

The shock modification can only be produced by the 
backreaction of the proton CR component. Therefore one 
can consider the extremely steep CasA spectrum in the 
radio range as indirect evidence of efficient nuclear CR 
production. The other important physical factor is the 
high magnetic field value which leads to substantial syn- 
chrotron losses of electrons with energies e c > 30 GeV: 
a magnetic field strength B s h = 200/xG is neded to re- 
produce the observed radio and X-ray synchrotron fluxes. 
Since we assume that the postshock field B 2 = <iB s , the 
downstream magnetic field amounts to about « £? 2 ~ 
1 mG in the shell (Fig.Uti) . This is roughly co nsistent with 
previous estimates (e.g. lAtovan et al.ll2000b|) . 

Such a large magnetic field B (r) considerably exceeds 
the value for a red supergiant wind, assuming that its typ- 
ical surface field at a radius of about 3 x 10 13 cm is about 
1 G. We believe that the large field strength in the shell 
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is likely either due to turbulent amplification of the red 
supergiant wind field by the shell formation in the final 
Wolf-Rayet phase or due to conside rable field amplifica - 
tion near the shock by CR streaming l|Lucek fe Belj200o|) . 
If the second factor is relevant, t he expected amplifie d field 
goes like B (R S ) oc ^N g (R s ) (jBell fc Lucekll200l|) . Our 
field Bq(R s ) is distributed according to this relation. The 
deduced field amplification is still well within the upper 
bound se t by the overall e nergy and momentum balance 
relations l|Volk et alJl2002t) . 

Our results show that even at the current epoch (when 
the SN shock propagates through the free red supergiant 
wind) the observed synchrotron emission is still deter- 
mined by the electrons accelerated during the shock prop- 
agation through the shell. Therefore the magnetic field 
B${r) in the wind zone is not a very relevant parameter 
for the fit to the observations. 

According to Eq. 10), we note also that the exponen- 
tial cutoff of the overall synchrotron spectrum at y max ~ 
10 18 Hz is roughl y independent of th e magn etic field, as 
pointed out bv lAharonian fc Atovanl (^999) for the case 
of a uniform field. Indeed the emission with frequency 
v ~ t'max is produced by electrons with momenta p ~ 
from their cutoff region. Since the maximum power of the 
synchrotron emission of electrons with momentum p is 
emitt ed at frequency v oc ctbBqP 2 (e.g. iBerezinskii et all 
1199(1 . and taking into account that p° nax oc V s / {a B ^B~ Q ) 
we conclude that v max oc V^/ctb, independent of B. A 
roughly 20% increase in f ma x would be required to bring 
the overall X-ray data into close agreement with the con- 
temporaneous model calculations (dashed curve), while 
leaving the low-frequency ranges unchanged at all cpoches 
(the low-energy electron spectrum is independent of the 
cutoff energy). Whereas the spectral hardening beyond the 
radio range is a basic feature of nonlinear shock acceler- 
ation, this is not true for the precise value of the overall 
high-energy cutoff. This cutoff depends on the detailed 
characteristics of the scattering wave field at long wave- 
lengths. With the above approximation of the diffusion 
coefficient by the Bohm limit, we see that a 20% increase 
in the characteristic value of V^/as is well within the 
variations of these quantities over the shell region (Fig[3J). 
Thus we can consider the agreement between theoretical 
results and observations as remarkably good. 

3.4. The secular decline of the radio flux 

Important information about the circumstellar parame- 
ters comes from the measured secular decline of the radio 
flux which is presented in Fig. 21 In the case of a very 
young SNR the population of CRs would be expected to 
be an increasing function of time, if the remnant was ex- 
panding into a homogeneous medium with constant den- 
sity and magnetic field. For CasA, even the existence of 
the decline of the radio flux itself is a strong indication 
that the SN shock at the current evolutionary stage prop- 
agates through a region of monotonically decreasing mag- 



7 




8 10 12 14 16 18 20 



log v, Hz 

Fig. 4. The rate of secular decline of the total flux of syn- 
ch rotron radiation. Th e diam ond shows the recent result 
of lO'Sullivan fc GreeiJ lll999l) at 15 GHz. All other data 
points are from lReej (J1990) . 

netic field strength B (Rs) and density Pn(R«). The mos t 
natural assumption, proposed bv lBorkowski et afl lll996h . 
identifies this region with the wind region. If we were to 
suppose that the CR population and the radio emission 
from its electron component are dominated by particles 
accelerated in this free wind region, the synchrotron flux 
would scale as 

S v oc £ (i? s )^ +1 '/ 2 AVi? s ), (14) 

where we represented the spectrum of radio emitting elec- 
trons in power-law form 

N c (p) = N (p/ Po y\ (15) 

In the wind region we have Bq(R s ) oc R~ , the number 
of accelerated electrons scales as the amount of swept-up 
mass Nq (x R s , and therefore 

S„ oc ^-(t- 1 )/ 2 . (16) 

In our case V s s» 2400 km/s, R s « 2pc and 7 w 2.6. For 
the secular flux decline 

dt 2 R s { ' 

we would therefore have dlgS v /dt ~ — 0.1%/yr. Its abso- 
lute value is considerably smaller than the observed one 
(Fig.©. 

The explanation for the large observed value of the 
secular decline of the radio flux comes from the idea that 
the radio emission is still dominated by the contribution 
of electrons accelerated at the previous stage during shock 
propagation through the dense shell. In this case the elec- 
trons undergo adiabatic cooling because of the expansion 
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of the downstream region, and one can rewrite relation 
(14) in the form 



S v oc B^'\-^\ 



(18) 



where a = (pi/p) 1 ^ 3 is the adiabatic factor which describes 
the particle adiabatic cooling, p is the gas density in the 
region where electrons are confined, and p\ is its value 
at the epoch when these electrons were accelerated. It is 
natural to assume that the downstream magnetic field is 
frozen into the gas. In this case it changes according to 
the relation B& oc a~ 2 . Let us assume that at the cur- 
rent epoch every downstream volume element has a oc R s . 
Then the expected decline of the synchrotron flux is 



dt 



-(27 + 3) 



(19) 



which gives dlgS v /dt f» — 0.8%/yr. This value is in the 
observed range dlgS^/dt w — 1.0%/yr to — 0.6%/yr (see 
Fig.0J). We take also into account that a revision of the 
data made be Rees (1990) led him to the conclusion that 
the secular decrease of the flux of CasA is about —0.8%/yr 
and is independen t of frequency. This is con sistent with 
the recent result of lO' Sullivan fc Greenl l|l999ft which gives 
—0.6%/yr at v = 15 GHz. We therefore conclude that the 
observed relatively large secular decreas e of the flux from 
Cas A is consistent with the assumption llShklovskvlll968t) 
that the flux is dominated by the contribution of electrons, 
produced in a previous epoch, which at the current stage 
undergo adiabatic cooling. Together with the decrease of 
the downstream magnetic field this cooling causes the rel- 
atively large rate of decrease of the radio flux. 

The picture described is reproduced in our model 
(Fig.QJ). One can see that the calculated secular decrease 
of the flux at v^lO GHz is about — 0.5%/yr and thus very 
close to the observations. At larger frequencies v> 10 3 GHz 
the theoretically calculated decrease has a rather compli- 
cated dependence upon v and is on average larger than 
at low frequencies. The highest secular decline predicted 
at the largest frequencies v> 10 16 Hz is due to synchrotron 
losses which lead to the decrease of the electron maximum 
energy. 

Since the shock currently propagates through the free 
red supergiant wind, the contribution of shell electrons to 
the total synchrotron flux will decrease in time compared 
with that of the electrons produced by the shock in the 
free wind. Therefore the rate of secular decline of the total 
synchrotron flux is expected to decrease from the current 
value —0.5%/yr towards — 0.1%/yr. 

3.5. The spatial structure of CasA 

The dowstream spatial distribution of radio emitting elec- 
trons is qualitatively consistent with the observed struc- 
ture of CasA. The brightest part of the remnant in our 
model is the swept-up shell. It is currently in the down- 
stream region near the contact discontinuity R p which 
separates ejecta and swept-up matter. We believe that it 



corresponds to the observed bright radio ring. The 'dif- 
fuse plateau' then corresponds to emission coming from 
the swept-up free wind matter further out. 

The small scale irregularities of t he radio emission, the 
so-ca lled compact radio knots (e.g. lAndersen fc Rudnickl 
Il99fih discussed earlier, which can not be reproduced in 
our spherically symmetric model, could well be due to 
amplification of the magnetic field at the bow shocks 
driven ahead of the so-called fast moving knots, in the 
form of dense clumps of SNR ejecta. Electrons with en- 
ergies e c < 1 GeV, which produce the radio emission, are 
strongly connected with the moving gas even on such small 
scales because of their small diffusion coefficient. Therefore 
their concentration also increases at the bow shock like 
the gas density. Together with the field amplification this 
leads to a strong increase of the radio emission from these 
relatively small volumes. 

We note that the field amplification leads to a de- 
creasing effective energy e e oc y/ v j B of the electrons 
which produce the synchrotron emission at a given fre- 
quency v. Together with the adiabatic electron heating, 
this leads to a brightness increase, and in parallel to a 
steepening of the radio spectrum because of the concave 
shape of electron spectrum in the corresponding energy 
range. Such a type of correlation between the knot bright- 
ness and synchrotron spect ral index is indeed observed 
l|Andersen fc Rudnicklfl996|) . 

The effect is expected to be larger in the outer down- 
stream region occupied by the swept-up free wind mat- 
ter and recently accelerated electrons compared with the 
swept-up shell region, where electrons after their produc- 
tion have already adiabatically cooled by expansion of the 
medium. Therefore a steeper synchrotron spectrum is ex- 
pected from the knots situated in the outer diffuse ra- 
dio plateau surrounding the bright ring, than from the 
knots in the ring. This is also in dicated by the observa- 
tions l| Andersen fc Rudnicklll996]) . 

3.6. High energy 7 -ray emission from CasA 

FigureEJrepresents the expected integral 7-ray energy flux 
components from NB, IC scattering on the background ra- 
diation field (cosmic microwave + optical/infrared), and 
hadronic collisions of CR protons with gas nuclei, respec- 
tively. 

As already mentioned, in contrast with our spherically- 
symmetric model one has to expect that not every part 
of the shock surface efficiently injects and accelerates 
CRs. At the local portions of the shock, which are cur- 
rently quasiperpendicular the suprathermal particle in- 
jection is presumably suppressed. Therefore the number 
of accelerated protons calculated within the spherical ap- 
proach should be corrected by s ome renormalization f actor 
/ rc < 1. It was argued earlier l|Berezhko et al-lfeoO^ that 
values / re = 0.15 to 0.25 are consistent with theoretical 
considerations (Volk et al., in preparation), the average 
requirements of the Galactic CR energy budget, and the 
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Fig. 5. Inverse Compton (IC, dash- dotted), nonthermal 
Bremsstrahlung (NB, dashed), and 7r°-decay (solid) inte- 
gral 7-ray energy fluxes as a function of 7-ray energy. The 
dotted line represents the same 7r°-decay 7-ray flux with 
an exponential cutoff exp(— e 7 /e™ ax ), wh ere e" ax = 4 TeV. 
The 1 TeV data point is from HEGRA l|Aharonian et all 
200 1]), the EGRET u pper limit for e 7 > 100 MeV is from 
Esposito et ail <|l996h . 



observed structure of the remnant SN 1006. One of these 
arguments is that for a typical SNR our spherically sym- 
metric model predicts a significantly larger CR produc- 
tion E c — 0.3 to 0.6-E sn than req uired for the Galactic 
energy budget Ej. ~ 0.1.E sn (e.g. iBerezhko et al 1 ll996t 
IBerezhko fc VoUdll997l boOOb)) . The 7r°-decay 7-ray flux 
presented in Fig. was calculated with a renormalization 
factor / rc = 1/6. Note that the number of accelerated elec- 
trons, which is required for the observed synchrotron flux, 
remains the same. Therefore the actual electron to proton 
ratio amounts to K cp = 0.02 instead of 4 x 10~ 3 . This also 
means that CRs absorb only E c — 10 49 erg at the current 
epoch. 

As one can see from Fig. [SI at e 7 = 1 TeV the 7r°-decay 
energy flux has a value of e 7 F 7 ~ leVcm _2 s _1 , which 
is close to the e xperimental value of th e reported integral 
TeV 7-ray flux l(Aharonian et al.ll200lil . Also the slope of 
the 7r°-spectrum lies in the reported range of allowed en- 
ergy spectra, which is quite large due to statistical uncer- 
tainties. The 7r°-decay energy flux hardens slightly above 
several GeV; nevertheless its value e 7 -F 7 remains practi- 
cally constant between 300 MeV and several TeV, and has 
a cutoff energy of about 30 TeV. 

On the other hand, the IC and NB fluxes at TeV 
energies are only at levels of e 7 F 7 ~ 10 _2 eVcm _2 s _1 
and below. Note that the electron scattering off the in- 
frared/optical background contributes about 50% to the 
IC 7-ray flux at energies e 7 < 1 Te V. This is by a facto r 
of two larger than the estimate of iGaisser et alJ lfl9 98). 
due to the fact that our electron spectrum is essentially 
steeper; for such a spectrum, the contribution of back- 
ground photons with higher energies becomes more im- 
portant. 



The cutoff energy of the IC spectrum is about 2 TeV. 
At 30 TeV, the leptonic 7-ray emission is entirely negli- 
gible. The contribution from the CMB target field drops 
because of the cutoff of the electron spectrum; the scatter- 
ing off the IR/optical photons is suppressed as a result of 
the decline of the Klein-Nishina cross section. The cutoff 
energy of the IC spectrum resulting from the CMB pho- 
tons is a bit lower than in the case of IR/optical radiation. 
Therefore at e 7 > ITeV the contribution of IR/optical 
background becomes progressively larger compared to the 
CMB photons, so that at e 7 = 10 TeV it exceeds the CMB 
contribution by a factor of 200. 

Figure [S] also shows the integ ral upper limit aboy e 
100 MeV as reported by EGRET l|Esposito et al.lll996h . 
This limit is well above the 7r°-spectrum; even an energy- 
resolved analysis of the EGRET data - which is not avail- 
able in the literature yet - is not expected to yield upper 
limits which would violate the model prediction. 

Pion production by the accelerated nuclei leads not 
only to 7-rays but also to secondary electrons and 
positrons. However, even for the large gas density 
N g < 50 cm~ 3 in the shocked shell, it can be easily shown 
that the ratio of secondary electrons to protons is more 
than two orders of magnitude lower than the value K cp = 
0.02 for primaries deduced here. Therefore the contri- 
bution of secondaries to the synchrotron and Inverse 
Compton emission and the nonthermal Bremsstrahlung 
is negligible. 

3.7. Remaining model uncertainties 

We believe that there is no significant difference between 
the calculated and measured TeV 7-ray flux of CasA. 
Nevertheless we have to ask ourselves to which extent our 
calculated 7-ray flux is a robust result and which effects 
could modify it. 

This concerns first of all the magnetic field configu- 
ration in the shell and in the free wind region. The RSG 
wind is most probably driven by strong wave activity from 
the star's outer convection laye rs, essentially Alfven waves 
llHartmarm McGregorllTflSrl. 

Their wavelengths should be quite large and unrelated 
to the gyro radius scale of the accelerated particles. The 
SNR shock propagates through this turbulent wind and 
injection occurs at the quasi-parallel portions of the mag- 
netic field lines. For large enough wave amplitudes in the 
wind, the mean field direction is no more relevant for 
the acceleration in a strong shock that produces in addi- 
tion its own large-am plitude resonantly scattering waves 
l|Lucek fc Belli l200f)fl . In the shell the situation is even 
more pronounced because of the compression by the Wolf- 
Rayet wind from behind and the result ing hydrodynamic 
instabilities ijGarcia-Segura et al.lll996|) . We therefore be- 
lieve that it is justified to use the Bohm limit for the 
diffusion coefficient, renormalizing the flux since the in- 
jection at the instantaneously quasi-perpendicular shock 
positions is reduced. Thereby we neglect the possible sys- 



E. G. Berezhko et al.: Gamma-rays from CasA produced by accelerated CRs 



10 

tematic effect of the mean field that in principle has a 
Parker spiral type configuration with the tendency to fur- 
ther reduce the injection rate. 

The second physical effect which can play a role in the 
current evolutionary phase is CR escape during the more 
recent phase of the SNR evolution. Let us assume that the 
red supergiant wind speed V w « lOkm/s does not exceed 
the Alfven speed, and that the field strength in the free 
red supergiant wind B W 2 is as low as 10 fiG. In this case 
the maximum momentum of accelerated protons 

R s V s /K(mc) (20) 

is at least 20 times smaller compared with its value 
in the shell since the minimum CR diffusion coefficient 
k cx 1/B increases by a factor of 20. In such a situation 
one should expect that the diffusion coefficient of particles 
with p > p™ ax , where p™ ax is the maximum momentum of 
CRs accelerated in the red supergiant wind region r > i?2 , 
exceeds the Bohm limit because the shock becomes unable 
to accelerate them and they consequently do not produce 
high level Alfvenic turbulence near the shock. It is like if 
the shock had become " old" , being unable to confine the 
highest energy particles it accelerated previously in the 
shell. Therefore these particles leave the SNR very soon 
after the SN shock reaches the red supergiant wind re- 
gion, resulting in a much smaller production of 7r°-decay 
7-rays with the highest energies compared with the previ- 
ous cases. For the electron X-ray emission this is a small 
effect since their main synchrotron radiation comes from 
the shell; their IC emission remains essentially the same 
with or without escape. 

To illustrate this scenario we present in Fig.El a 7- 
ray spectrum which compared to the previous one has an 
additional cutoff factor exp(— e 7 /e™ ax ); the cutoff energy 
£ max wag ge |- |- Q 4 rp e y w j 1 j c j 1 corresponds to the escape of 

particles with energies above 50 TcV, lowering the proton 
maximum energy by a factor of 20. As one can see from 
Fig.[5]this would still be consistent with the HEGRA flux, 
but steepen the spectrum at 1 TeV. 

The integral NB and IC energy fluxes are almost 2 or- 
ders of magnitude lower than the corresponding 7r°-decay 
flux. Even taking into account the local thermal far in- 
frared fluxes can not increase the IC flux by more than a 
factor of 2 (R.J. Tuffs, private communication). Therefore 
the dominance of the 7r°-decay flux is a rather robust re- 
sult. 

4. Summary 

Perhaps the most important result of our considerations is 
that the spectral shape of shock accelerated electrons with 
their essential synchrotron cooling in the downstream re- 
gion is very well consistent with the observed synchrotron 
emission. To reproduce a very steep radio spectrum S v oc 
v~ - 77 the shock must be strongly modified. This shock 
modification can only be produced by accelerated protons 
if they are also efficiently injected into the acceleration 
process, as it was assumed in the calculation. 



The significant synchrotron losses of electrons in the 
strong interior magnetic field makes their spectrum steep 
N c cx e~ 3 also at high energies e c > lOGeV. This leads 
to a flat connection of the spectral energy distributions of 
the observed radio and X-ray synchrotron emissions. 

The rather high secular decline of the synchrotron radi- 
ation observed in the radio range is naturally reproduced 
in our model, because this range is dominated by elec- 
trons which were accelerated at a previous epoch and un- 
dergo currently adiabatic and synchrotron cooling in the 
expanding downstream region. 

We find that after reduction of the predictions of the 
nonlinear spherically-symmetric model by a renormaliza- 
tion of the number of accelerated nuclear CRs, to take 
account of the large areas of quasiperpendicular shock re- 
gions of a SNR, good consistency with all observational 
data can be achieved, inclu ding the reported TeV 7-ray 
flux l|Aharonian et alJ l200lf) . The used renormalisation 
factor is consistent with the need that the average ac- 
celeration efficiency of a typical SNR within our model 
scenario should meet the requirements for Galactic CR 
acceleration. 

In addition, our calculations show that at all energies 
above 1 GeV the 7-ray production is dominated by ir°- 
decay. At TeV energies the expected 7r°-decay flux ex- 
ceeds the IC and NB fluxes by a factor of about seventy. 
Therefore the leptonic emission is totally inadequate to 
explain the observed TeV 7-ray flux. The 7r°-decay spec- 
trum oc e" 1 extends up to 30 TeV, whereas the IC and 
NB 7-ray fluxes have a cutoff at about 1 TeV. Therefore 
the detection of 7-ray emission at 10 TeV and above would 
imply further evidence for its hadronic origin. 

We conclude that the observed properties of the radio 
and X-ray emission can be explained within the assump- 
tion that the SN blast wave is the main source of energetic 
particles in Cas A. The CR production efficiency and the 
electron to proton ratio implied by these multi-wavelength 
observations are consistent with the requirements of the 
nuclear CR sources in the Galaxy. 
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